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Abstract. The trade-off between the information gain and the state disturbance is derived 
for quantum operations on a single qubit prepared in a uniformly distributed pure state. The 
derivation is valid for a class of measures quantifying the state disturbance and the information 
gain which satisfy certain invariance conditions. This class includes in particular the Shannon 
entropy versus the operation fidelity. The central role in the derivation is played by efficient 
quantum operations, which leave the system in a pure output state for any measurement outcome. 
It is pointed out that the optimality of efficient quantum operations among those inducing a given 
operator-valued measure is related to Davies' characterization of convex invariant functions on 
hermitian operators. 



1. Introduction 

Disturbance caused by an attempt to gain classical information about the quan- 
tum state of a physical system is one of the essential features of quantum theory 
[I1G1I31G1E1E1IIHEI The balance between the information gain and the state dis- 
turbance can be quantified in many different ways, depending on a specific physical 
scenario. Analytical derivations of the exact form of such trade-offs are usually a 
non-trivial matter, as the optimization needs to be carried out over all possible 
quantum operations that can be applied to the system. Nevertheless, in certain 
instances, for example the case of estimation and operation fidelities for a single 
d-level system in a uniformly distributed pure state, optimization by analytical 
means has been shown to be possible 0. The derivation of this last result was 
enabled essentially by the fact that both the operation and the estimation fidelities 
are quadratic in the operators that constitute the Stinespring decomposition of a 
quantum operation. 

The situation becomes considerably more involved when the parameters quan- 
tifying the state disturbance or the information gain have a more complicated 
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dependence on the quantum operation. A good example is the gain of information 
expressed in terms of the Shannon entropy. This quantity arises naturally in the 
framework of Bayesian inference, when we try to quantify how well the a posteriori 
probability distribution "pins down" the measured state (H1E3^3- In this paper, 
we present a derivation of the trade-off between the information gain and state dis- 
turbance, when the physical system under consideration consists of a single qubit 
prepared in a pure state, and both the information gain and the state disturbance 
are quantified using parameters satisfying certain invariance criteria. As a concrete 
example, we will consider the trade-off between the Shannon entropy as a measure 
of the information gain combined with the state disturbance quantified in terms of 
the average operation fidelity. 

An intermediate step in this derivation requires demonstrating that it is suffi- 
cient to consider the so-called efficient j^] (also known as Luders-type f4| or ideal 
|12j ) operations, for which obtaining a specific measurement outcome leaves the 
system in a pure state. Physically speaking, efficient quantum operations should 
introduce minimum disturbance while inducing a fixed generalized measurement, or 
equivalently, maximize the amount of classical information for a given Stinespring 
decomposition of a quantum operation. While the latter property is straightfor- 
ward to prove for the Shannon information gain or the average estimation fidelity, 
the former property in the case of the operation fidelity has been conjectured a 
while ago by Barnum 0], and proven by him only recently jH] using the convex- 
ity of certain operator maps involving tensor products |13j . We show here that 
in the case of uniform distributions generating unitarily invariant measures of the 
state disturbance, the optimality of efficient operations can be analyzed following 
an alternative route, based on results obtained in the field of convex analysis on 
operators ^3J DI3 EI This connection is valid for arbitrary finite- dimensional 
systems, and it is an interesting question whether any of these results can be carried 
over to Hilbert spaces of inifinite dimension. 

This paper is structured as follows. First we discuss the convexity properties 
of measures characterizing information gain and state disturbance in Sect. 2 and 3 
using the examples of the Shannon information gain and the operation fidelity. In 
Sect. 3 we also link the optimality of efficient quantum operations to the convexity 
of unitarily invariant functions defined on hermitian operators. In Sect. 4 we derive 
the trade-off between the information gain and the state disturbance for a single 
qubit. Sect. 5 gives two other measures of information gain and state disturbance 
which fall under the general formalism used in the preceding section. Finally, Sect. 
6 concludes the paper. 

Before we pass on to detailed discussion, let us first summarize briefly the ele- 
mentary facts needed further in this paper. The subject of our interest is a quan- 
tum operation, i.e. a trace-preserving completely positive map acting on density 
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operators g. The map has the Stinespring decomposition of the form 

Q^Ys^QAU- (!) 



with the family of the operators {A ril } satisfying the condition 

E44 = i (2) 

r/j, 

which guarantees the preservation of the trace of the density operator. We assume 
that the index r describes the classical outcome of the measurement, while the 
summation over the index [i is responsible for additional averaging due to the 
imperfections of the measuring apparatus, resulting in a decrease of classically 
available information. If we are interested only in the classical outcome of the 
measurement, then it is sufficient to consider a positive operator- valued measure 
{M r } induced by the quantum operation {A ril } according to 

^ = E^V (3) 

The probability of obtaining a specific result r is given by Tr(M r g), and the out- 
come r is associated with a conditional transformation of the density matrix given 

by 

'"^ki)^ Al - (4) 

The gain of classical information depends only on the operator-valued measure 
{M r }, whereas the state disturbance is affected by the specific form of the quantum 
operation. Of course, there are many quantum operations corresponding to a given 
operator-valued measure {M r }. Given an operation {^4^} we can always augment 
the associated information gain by assuming that both the indices r and \x are 
available classically. Conversely, we can ask which quantum operation inducing the 
generalized measurement {M r } minimizes the state disturbance. A good candidate 
for this operation is {\/ M r }, as it retains the purity of the input state for a given 
measurement outcome, and also preserves the relative phases of the input state in 
the basis diagonalizing M r . As it will be discussed in Sect. 3, this intuition turns 
out to be correct when the measure of the state disturbance can be represented as 
a sum of contributions given by the values of a certain function on the operators 
A-t^Ar/j,, when this function is unitarily invariant and concave on the subspace of 
diagonal operators. 
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2. Shannon information gain 

As an example of a measure of classical information gained from a quantum 
operation, we shall consider here the decrease of the Shannon entropy between the 
a priori and a posteriori probability distributions defined on the manifold of pure 
states. Let us denote by dip the normalized invariant integration measure on the 
manifold of pure states. This measure is induced by the Haar measure associated 
with the Lie group of unitary transformations on the corresponding Hilbert space. 
We will assume that the initial a priori probability distribution is given with respect 
to this measure by a function p(ip). Therefore the initial entropy is given by the 
expression 

H\p{ip)\ = - f ctyp(V01og 2 pM- (5) 

Obtaining a specific outcome r of the measurement yields a conditional a posteriori 
probability distribution p{ip\r) calculated according to the Bayes' rule as 

p(r\ip)p(ip) 

p^V) = —rff) — (6) 

where p{r\ip) = {ip\M r \ip) is the standard quantum mechanical probability of ob- 
taining the outcome r for an input state \ip), and p(r) is the average: 

p(r) = J dip p{r\ip)p{ip) (7) 

The Shannon entropy of the conditional a posteriori distribution p(ip\r) is given 
by the expression H\p(ip\r)\. As we are interested in the average decrease of the 
entropy after the quantum operation is carried out, we need to take the difference 
H\p(ip)] — H\p(ip\r)] and average it over all possible outcomes of the experiment 
with the probability distribution p(r). Thus the final expression for the average 
gain of the Shannon information H from a measurement described by a positive 
operator-valued measure {M r } is given by: 

H = ^p(r){i7[p(V0]-#b(V'k)]} 

r 

= E/ d'/'P(''l«)p(</')log 2 ^ i 

In the following, it will be convenient to represent the average information gain as 
a sum 

H = Y,n{M r ), (9) 
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where H(M) is a function defined for semipositive operators M according to: 

H(M) = I d0p(VOMJft|V)log a r,/! 1 !!'!, , i (10) 

Let us note that H(M) is positively homogeneous of degree one. 

Physically, we anticipate that combining two elements M\ and M 2 of the oper- 
ator measure into one operator M\ + M 2 should result in the loss of information. 
Such a loss of information should correspond to the inequality 

H{M X ) + H(M 2 ) > H(Mi + M 2 ), (11) 

which is equivalent to the convexity of H(M) due to the fact that TL(M) positively 
homogeneous of degree one. The above inequality is a general property of mutual 
information [18J, and it follows elementarily from the convexity of the function 
x log 2 x on the interval [0, 1] : 

txi log 2 x\ + (1 — t)x2 log 2 x 2 > [txi + (1 - t)x 2 ] log 2 [txi + (1 - t)x 2 ] (12) 

Inserting 

t _ P(n) ^ = p{ri\ip) ^ = p{r 2 \4>) - 13 . 

p(n) +p(r 2 y p(ri) ' p(r 2 ) 

yields the convexity of the integrand in lfTT)|) . and integrating it with j d'0p(t/') 
proves the convexity of H. By mathematical induction, we therefore obtain that 
for a quantum operation {^4 r ^} the information gain is maximized by assuming 
that both the indices rfj, are available classically, and by constructing a "finer" 
generalized measurement from the quantum operation as M r/1 = At^Ar^. An 
analogous statement holds also for any other measure of information gain which is 
decomposable to the form Q with a positively homogeneous convex function TL. 



3. Operation fidelity 

If the initial state \tp) of the system is chosen according to a probability distri- 
bution p(V0) then the operation fidelity F for a quantum operation {A rfl }, defined 
as the average projection of the final density matrix onto the initial pure state, is 
given by the expression: 

F = Y1 [wpWMArnW 2 . (14) 

For a uniform distribution of input states with = 1 the integral f dtp can 

be carried out analytically in a d-dimensional Hilbert space (U EH], yielding the 
following expression for the mean operation fidelity [H[5]: 

F = d(dTl) D^/A*) + Mr/]- (15) 
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By applying the singular value decomposition to operators A r ^, it can be shown 
by elementary means that TrA r ^ < Tta/ A|„A r u, and consequently the operation 



fidelity can only increase when replacing A rfl with \J At^Ar^, while the general- 
ized measurement induced by the operation (and consequently any measure of the 
information gain) remain the same. This implies that it is sufficient to restrict our 
considerations to quantum operations composed of semipositive definite hermitian 
operators A rfJ- . In this case, we can represent the operation fidelity as a sum over 
rfi of a certain function T of the products At^A r/1 : 

F = Y,HAl,A r ,). (16) 

The function T(M), given explicitly by: 

HM) = -^r^r (TrM + \Tv\/^\ 2 ) (17) 

is well defined for all semipositive hermitian operators M, which form a convex 
set. Let us note that similarly to the function TC(M) introduced in the preceding 
section it is also positively homogeneous of degree one. 

We would like now to find a quantum operation which maximizes the mean op- 
eration fidelity under the constraint of inducing a specified operator- valued measure 
according to (J3J|. Recalling the introductory discussion, we expect that the optimal 
operation is given by square roots of the operator measure {y/ M r }. In order to ob- 
tain this result, it is sufficient to show that for any semipositive definite operators 
Mi and Mi the function T{M) satisfies the following inequality: 

JF(Mi) + T{M 2 ) < T(M X + M 2 ), (18) 

which again is equivalent to its concavity, as T{M) is positively homogeneous. 
Assuming that lfT%|l holds, we would immediately obtain that for any two elements 
A rfJ-1 and A rfl2 of a quantum operation we have: ^(At^Aj.^) + T{A} ril2 A rpL2 ) < 
^(At^Ar^ + Ai^ 2 A r/12 ), and by induction: 

Y^H^A^) < T [S^A^ArA = F(M r ), (19) 

M \ (U / 

which would prove our thesis. 

In order to demonstrate lfT8|) . let us first note that T is invariant with respect 
to unitary transformations of its argument, i.e. T{W Mil) = ^{M) for any semi- 
positive M and unitary U . Therefore ^(M) can depend only on the eigenvalues 
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of M. Let us consider a function <fr denned on d- dimensional vectors u G with 
nonnegative coordinates, given by the value of T on a diagonal operator diag(u): 

0(u) = ^(diag(u)). (20) 

In the case of T given by (|TK|l . the function (ft reads: 

M = ^ TT) (E^+(E^) ) (21) 



d(d- 



i=0 \i=0 



and it is symmetric, i.e. invariant with respect to the permutations of the coordi- 
nates Uo,ui, . . . , Ud-i of the vector u. It is easy to verify that for any pair u, v the 
following inequality holds: 

0(u) + 0(v) <0(u + v), (22) 

which analogously as before means that (ft is concave due to its positive homogene- 
ity. Indeed, by writing explicitly the left and the right hand sides of the above 
inequality we have: 




J \ui + Vi)(uj + Vj) (23) 



and it is straightforward to see that the inequality holds separately for every term 
with fixed % and j, which proves 1)22)1. The inequality l(2*2*|) immediately implies the 
inequality (|T8|) for pairs of operators M\ and M2 which commute and therefore can 
be diagonalized simultaneously in the same orthonormal basis. 

The definition l)20)l introduces a correspondence between a unitarily invariant 
function T and a function (ft which is symmetric, i.e. invariant with respect to the 
permutations of the coordinates of its argument. The crucial step of the reasoning is 
the application of a theorem proven first by Davies [2] and by Friedland ^] stating 
that in this setting the convexity (concavity) of T is equivalent to the convexity 
(concavity) of cj> on their respective domains. Hence (|2*2"|) implies the inequality (fTH|) 
also for non-commuting pairs of operators M\ and Mi. This immediately proves 
that for a given generalized measurement {M} the quantum operation minimizing 

the state disturbance quantified with the mean operation fidelity is given by {\f~M~} 
in the case of a uniform distribution on pure input states. 



4. Trade-off for a single qubit 



The results described in the preceding sections allow us to restrict the search 
for the information gain versus state disturbance trade-off to efficient quantum 
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operations generated from positive operator-valued measures {M r } by taking the 
square roots of its elements {\/ M r }. Furthermore, we will consider T and H 
which are respectively concave and convex unitarily invariant functions with the 
property of positive homogeneity. For a single qubit, these assumptions severely 
restrict the number of parameters characterizing the elements M r of the operator- 
valued measure which play a non-trivial role in the derivation. In fact, as we will 
see in a moment, the only relevant parameter is the ratio of the eigenvalues of M r . 
Let us start by using the positive homogeneity of T and TL to write: 

F(M r ) = iTr(M r )^ . 

2 v r; \TV(M r 



H(M r ) = -Tr(M r )TL 




(24) 

\^Tr(M r 



and denote £ r = Tr(M r )/2. Of course £ r > 0, and taking the trace of ^2 r M r = 1 
yields a summation condition on {£ r }: 

2> = i. ( 25 ) 

r 

Next, let us denote the eigenvalues of the renormalized operators 2M r /Tr(M r ) as: 



Tr(M r ) 



The positivity of the eigenvalues requires that — 1 < x r < 1, and if we further as- 
sume with no loss of generality that the eigenvalues are arranged in a non-increasing 
order, we can restrict our interest to the range < x r < 1. In the new variables 
we can now write: 

HMr) = Hrf(Xr) 

H(M r ) = irh{x r ) 

where the functions /, h : [0, 1] — > R are defined as: 

f(x) =.F(diag(l + x,l-x)) (28) 
h(x) =K(diag(l + x,l-x)) (29) 

It is straightforward to show that the concavity of T and the convexity of TL 
together with their positive homogeneity imply that the functions / and h are 
respectively concave and convex. Furthermore, the converse is also true. 

Physically, we expect that / is strictly decreasing, whereas h is strictly increas- 
ing on their domain [0, 1] . The reason for this is that the parameter x characterizes 
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the imbalance between the eigenvalues of M r for a fixed trace Tr(M r ), and bal- 
anced pairs of eigenvalues do not disturb the state, whereas unbalanced pairs of 
eigenvalues are responsible for the information gain. If we now assume that the 
function / is strictly monotonic, then the information gain H can be expressed as: 

H = = Y,tMrHfM)) = O r l ){f{Xr)) (30) 

r r r 

Let us now recall the notion of a concave envelope [2Hj of a real function % : C — > 
R defined on a convex set C. The concave envelope, which we will denote by 
X '■ C — ► ft is defined as a concave function such that x( x ) — x( x ) f° r an x £ C, 
and also that for any other concave function x '■ C — > R satisfying x{x) > x( x ) 
everywhere we have also x(x) > x( x ) everywhere on C. In the one-dimensional 
case the concave envelope is given explicitly by the expression 

X{x) = sup {tx{xi) + (1 - t)x{x 2 ) | < t < 1, txx + (1 - t)x 2 = x}. (31) 

t,Xl,X2 

In our derivation, we will take x = ho f _1 , and xi, X2 in the above formula are any 
two points from the image of /. The concave envelope x = h° f 1 can be used 
to estimate H in 1(30)1 from above, and furthermore the concavity of x allows us to 
apply Jensen's inequality. Combining these two steps yields: 



H < 

Thus finally 



£>X(/K)) < X (Y,trf(*r)j = X(F). (32) 



H<hof-\F), (33) 

where ho / _1 is defined according to fpTTf) . Let us note that in general the con- 
cavity of / and the convexity of h along with their monotonicity do not guarantee 
automatically that the composition ho f^ 1 is concave itself, hence the need to take 
its concave envelope. This can be seen most easily when both / and h are doubly 
differentiable on [0, 1]. Then we have: 



h" f - h' f" 

"'°r*Y = hf u * f (34) 



and the conditions h! > 0, h" > 0, /' < 0, and /" < do not imply a well-defined 
sign of the right-hand side of the above formula. 

The inequality 113311 can be easily saturated. If we take a two-element quantum 
operation of the form 



= diag(v/(l + z)/2, V(l-x)/2),i 2 (x) = diag(v/(l - x)/2, y/(l + x)/2)} 
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with < x < 1, we can generate any point of the graph of the function h o 
Furthermore, by taking suitable combinations of two quantum operations that have 
the above form 

{Vtii(xi), V~tA 2 (x l ),Vl-tA l (x2), Vl-tA 2 (x 2 )} 

with < t < 1 we can generate an arbitrary point of the graph of the function 
h o which follows from (pTTj). Let us also note that if the composition h o f^ 1 
is a strictly concave function, then the equality sign in (|3*3*|) is reached if and only 
if the parameters x r are equal. This means that for a strictly concave h o f^ 1 all 
the elements of an operator valued measure saturating the trade-off need to have 
the same ratio of their eigenvalues. 

Let us now specialize the above result to the trade-off expressed in terms of the 
operation fidelity versus the Shannon information gain. In the case of the average 
operation fidelity, the function / defined in (|2*%|) has the following explicit form: 

/(x) = i(2 + Vl-x 2 ) (35) 

Furthermore, as sketched in Appendix A, it is straightforward to obtain a closed 
expression for the Shannon information: 

/,(*) = — log 2 (— ) + - 2 log 2 (l -*>) - ^ (36, 

which is a particular case of the general result obtained by Jones [21]. In the case 
of / and h given by (|3*K|l and (|3l)|l the composition h o / _1 is strictly concave on 
the whole image of /, which is also shown in Appendix A. Consequently, inverting 
the function / given in 1|35|) brings the trade-off inequality: 

H < h{y/l- (3F-2) 2 ) (37) 

with h given by (|3T)|l . As the right-hand side is strictly concave over the domain of 
F, Jensen's inequality is saturated only for operations for which all x r are equal. 
This means that the ratio of the eigenvalues for all the elements of the operator 
measure {M r } needs to be constant across the index r. 

5. Other measures 

In order to illustrate the generality of the derivation presented in the preceding 
section, let us discuss two other measures of information gain and state disturbance 
for which the above reasoning leading to the trade-off for a single qubit holds as 
well. 

As an example of another measure of the information gain satisfying all the 
above properties let us recall the example of the mean estimation fidelity |22ll23ll2"3] 



[Author and title] 



11 



for which the trade-off against the operation fidelity has been studied in 
Using the notation of the present paper we can write the estimation fidelity for a 
uniform input distribution in d dimensions as a sum G = Y^ r QiMr) of terms given 

Q{M) = —L— [Tr(M) + || V^H 2 ] (38) 

where 1 1 • 1 1 stands for the standard Euclidean operator norm. The above function 
again is positively homogeneous and convex: 

g{M 1 )+g{M 2 )>Q{M 1 +M 2 ), (39) 

the latter property following immediately from the lower bound on the second term 
in l(3*H|l for a sum of two semipositive definite hermitian operators Mi and M 2 : 



Mi + M 2 || 2 = sup (0|Mi + M 2 \4>) = (<h\Mi + M 2 \<h) 
= (0o|Mi|0o) + (0o|M 2 |0o) 

40 

< sup {<j>\ Mi \<j>)+ sup (0|M 2 |0) 



'M 1 \\ 2 + \\dM 2 



2 



where \4>o) is the eigenvector of Mi + M 2 corresponding to its largest eigenvalue. 
Specialized to the qubit system, Q induces the function g: 

ff = g(diag(l + x,l-x)) = -(x + 3). (41) 

D 

As an alternative measure of the state disturbance we will consider the ab- 
solute value of the scalar product between the input and the output state, and 
its generalization to density matrices [23 EE! E3 In this paper we will call it 
the Bures-Uhlmann fidelity in order to distinguish it from the fidelity discussed 
in Sect. 3, which is a square of the former |28j . For a distribution on pure input 
states characterized by a probability distribution p(V0, the average Bures-Uhlmann 
fidelity is given by 

B = Y,t d^p(V)\/ W&^MMArMY (42) 

rfj, J 

The demonstration that quantum operations composed of semipositive definite her- 
mitian operators optimize the Bures-Uhlmann fidelity turns out to be substantially 
more complicated than in the case of F. The particular case of a single qubit with 
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p(ip) = 1 is discussed in Appendix B. If we restrict our attention to quantum oper- 
ations composed of hermitian and semipositive definite operators A r ^ 



we can write B as a sum: 



(43) 



r/j, 



where the function B is defined on hermitian semipositive operators by the equa- 
tion: 

B(M) = j d^p(iP)\J {^M^^^fil^). (44) 

The function B is positively homogeneous, and it is also concave for commuting 
operators M\ and M 2 . In order to show this, let us introduce an orthonormal 
basis \i) in which both the operators are diagonal: M = Y^t=o u i\^) {^\ an d ^2 = 
YdZlvi\i){i\, and denote pi = \{i\ip)\ 2 . By treating ^/ul, y^, JYllZlPjUj, and 



Yl'j=oPj v j as f° ur independent nonnegative numbers, it is straightforward to 



J 3 

show that 



d-i 



\ 3=0 



d-l 

^PjVj < yJUi + 
j=0 



d-l 



\ 3=0 



(45) 



This inequality, summed over i with the weights pi yields: 

(V'lv^lV')v / (V'|M 1 |^) + (^|v / M2|V')v / (V'|M2|V') 



< 



\\Jm 1 + M 2 \ij))\l (V>|M! + M 2 \ip) (46) 



which integrated over pure states with J dippfy) implies that for any two commut- 
ing operators Mi and M2 we have: 



B(Mi) + B{M 2 ) < B(M l + M 2 ). 



(47) 



Furthermore, if p{ip) = 1 then B is unitarily invariant, and we can again use 
Davies' result to show that B is concave for all, not necessarily commuting, pairs 
of semipositive hermitian operators. 

The above result means in particular that for a qubit prepared in a uniformly 
distributed pure state efficient operations composed of hermitian semipositive ele- 
ments are optimal from the point of view of the trade-off. In this case it is sufficient 
to consider a function b of a single real parameter x defined analogously to Sect. 4 
as: 



b(x) = £(diag(x + l,x - 1)) 
2 



15x 2 



[(1 + x 2 )Vl -x 2 + 7x 2 - 1]. 



(48) 
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F B 



Fig. 1: The quantum mechanical trade-offs between the information gain and the 
state disturbance with the state disturbance quantified as the average operation 
fidelity F and the average Bures-Uhlmann fidelity B, and the information gain 
measured using the change H in Shannon entropy, and the average estimation 
fidelity G. In all four cases the trade-offs are characterized by concave composite 
maps specified in the upper right corners of the graphs. 



The resulting trade-offs for any combination of F and B with H and G are depicted 
in Fig. ^ It is seen that in all four cases the composite maps characterizing the 
trade-offs are concave themselves and therefore equal to their concave envelopes. 

6. Conclusions 

We have discussed the convexity properties of measures quantifying the infor- 
mation gain and the state disturbance in quantum operations, using the examples 
of the average operation fidelity and the Shannon entropy. Such convexity prop- 
erties can be expected from all physically motivated measures, and in the case 
of uniform a priori distribution which implies invariance with respect to unitary 
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transformations, we can use the theory of convex invariant functions to analyze 
the convexity. In the case of quantum operations on a single qubit, the trade-offs 
can be described using a general inequality, resulting from the existence of a sin- 
gle relevant parameter characterizing the elements of the positive operator-valued 
measure through the ratio of their eigenvalues. 
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Appendix A 

In order to evaluate the function h defined in (|2H|l let us introduce the standard 
parameterization of the manifold of the single-qubit pure states: 

■ « = ^in«d S d„, (49) 

with < 9 < it and < tp < 2ir. In this parametrization we have: 
h(x) = H(diag(l + x, l-x)) 

1 r 

= — / d9 sin 9 (1 + x cos 9) log 2 (l + x cos 9) 

2 Jo 

" ^'^(^H 1 ^- 2 )-^. (50) 

In order to demonstrate the concavity of h o f~ x we need to check the concavity 
of the function h(\/l — x 2 ) on the interval x G [0,1]. The second derivative of 
/i(Vl — x 2 ) is given by 



— h(\/l -x 2 ) = — 

dx 2 ^ X) 4(l-X 2 )5/2 

(51) 

The expression in the square parentheses is equal to for x = 0, and then its 
derivative is equal to —4(1 — x 2 ) 3 ^ 2 /3x 2 In 2, which is strictly negative for x G]0, 1[, 
This implies that the expression under consideration takes negative values on that 
interval. Consequently d 2 h(\/l — x 2 )/dx 2 is negative for x e]0, 1]. 



/ 1 - Vl -x 2 \ (Ax 2 + 2)y/l -x 2 

lg2 ^rTTrT^f J + 3^2ln2 
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Appendix B 

Let us consider the contribution to the Bures-Uhlmann fidelity generated by 
an element of quantum operation of the form UA where U is unitary and A = 
diag(Vl + x, yfl — x) with x £ [0, 1] . It will be convenient to switch to the Bloch 
vector representation, where the density matrix of a pure state \ip){ip\ is represented 
by a three-dimensional real vector r according to: 

1 .(i + r r .*) (52) 



2 

where <x is a vector composed of the three Pauli matrices. If we denote A = 
diag(Vl — x 2 , \f \ — x 2 , 1) and a = (0, 0, x) T , then the expectation value (V>|^^4|V') 
can be written as: 

(tp\A^A\^) = 2(l + a T -r), (53) 

where the dot • denotes matrix multiplication, and the conditional transformation 
IV'}^! I— ► UA\ip)(ip\AW'/(ip\A 'A\ip) is given in the Bloch representation by 

O • ( A • r + a) 

r ^ — ^ ^— -L 54 

1 + a T • r y ' 

where O is the rotation of the Bloch vector corresponding to the unitary transfor- 
mation U . 

The integrand in lj%2*jl for p(ip) = 1 can now be written as 



(4>\AlA\iP}\(iP\uA\%b}\ = -^=J(l + aT • r ) 2 + (1 + a T • r)r r • O • (A • r + a). 

a/2 v 

(55) 

Following PHI, we will parametrize the Bloch vector r as 

r = (sin 9 cos <p, sin 9 sin (p, cos 9) T (56) 

and then the Bures-Uhlmann fidelity for the operator UA is the expression given 
in l|55p integrated over 

■ A = i (57) 

Let us first consider the second integral over the azimuthal angle (p, We can write 
the scalar product appearing in the integrand as: 

r T • O • (A • r + a) = Tr[0 • A(r + a) • r T ] (58) 

and use the Schwarz inequality for functions on [0, 2tt] to arrive at the upper bound: 

r 2w 



_t_ [ ZW dtp 

V2 



[ ^ — a/(1 + a T • r) 2 + (1 + a T • r)Tr[0 • A(r + a) • r T ] 
Jo 27r v 



1 



< V (1 + a T • r) 2 + (1 + & T • r)Tr(0 • A • V) (59) 
V2 v 
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where the matrix V is given by 

V= f ^(r + a) -r T = diag(± sin 2 ± sin 2 fl, x cos # + cos 2 #) (60) 
Jo ^ 

Let us now look for a rotation O that will maximize the right hand side of the 
inequality l[59|) . As O enters this expression only through the trace Tr(0 • A • V) 
and furthermore both A and V are invariant with respect to rotations about the 
z axis, it is sufficient to consider rotations of the form O = R z (a) • H x ((3). The 
trace written explicitly takes the form: 

Tr(0 • A • V) = - s/l - x 2 sin 2 9 cos a(l + cos (3) + (cos 2 9 + x cos 9) cos (3 (61) 

It is seen that this expression reaches the maximum for a = irrespectively of 
the values of any other parameters. Setting cos a = 1 and introducing a new 
integration variable t = cos 9 for the polar angle 9 we can write the expression for 
the Bures-Uhlmann fidelity in the form of an integral: 

~ f dt J 2(1 + xt) 2 + (1 + xt)[\/l - x 2 (l - t 2 ){l + cos (3) + 2{t 2 + xt) cos f3] 

(62) 

It can be verified by numerical means that for any value of x G [0, 1] this integral 
reaches its maximum value for (3 = 0, corresponding to U = 1. When f3 = the 
above integral can be evaluated analytically, with the final result given in 
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